Ultrahigh Energy Nuclei in the Turbulent Galactic Magnetic Field 



G. Giacinti ab , M. KachelrieB c , D. V. Semikoz a - d , G. Sigl b 

" AstroParticle and Cosmology (APC), Paris, France 
*//. Institutfiir Theoretische Physik, Universitdt Hamburg, Germany 
c Institutt for fysikk, NTNU, Trondheim, Norway 
" Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, Russia 



Abstract 

in this work we study how the turbulent component of the Galactic magnetic field (GMF) affects the propagation of ultrahigh energy 
heavy nuclei. We investigate first how the images of individual sources and of the supergalactic plane depend on the properties of 
the turbulent GMF. Then we present a quantitative study of the impact of the turbulent field on (de-) magnification of source fluxes, 
' due to magnetic lensing effects. We also show that it is impossible to explain the Pierre Auger data assuming that all ultrahigh 
ryy energy nuclei are coming from Cen A, even in the most favorable case of a strong, extended turbulent field in the Galactic halo. 
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1. Introduction 



(-h The composition of the cosmic ray (CR) flux above ener- 
Q_igies E > 10 17 eV is a topic of current debate. The shape of 
' the CR spectrum |Q1 0, O at the highest energies is compat- 
5_i ible with either a proton or a heavy nuclei composition. On 
the one hand, the correlation of ultrahigh energy cosmic rays 
C$ (UHECRs) arrival directions with the large scale distribution of 
active galactic nuclei reported by the Pierre Auger collabora- 
tion (4J, |5J could be consistent with a composition dominated 
by protons. On the other hand, the possibility to accelerate nu- 
clei to higher energies as protons may favour a transition to 
heavier nuclei at the end of the UHECR spectrum. Meanwhile, 
various results of measurements of the UHECR composition 
have been presented: The Pierre Auger Observatory air shower 
measurements suggest a gradual shift towards a heavier compo- 
sition starting at a few times 10 18 eV up to 4 x 10 19 eV above 
which the statistics is currently still insufficient for composition 
studies ||6|]. The Yakutsk EAS Array muon data are compati- 
ble with this result [7]. However, both HiRes measurements 18J 
and preliminary studies of the Telescope Array ii do not agree 
with this shift. The results of their composition studies favour 



> 



O 



proton primaries. 

In this work, we extend our previous analysis of the implica- 
tions of a possible heavy primary composition on the UHECRs 
propagation in the Galactic magnetic field (GMF) 01 Oil : We add 
a turbulent component to the regular GMF and investigate in 
details its influence. We present in this paper simulations for 
iron nuclei with the energy E — 60 EeV, since this is the high- 
est energy where present experiments collect still a reasonable 
number of events. 

T. Stanev suggested one of the first models of the regu- 
lar GMF including a spiral arm-like structure of the field in 



the Galactic disk 111 111 . Other models were later proposed 
in Refs. HI Q. The authors of Refs. [14] and [15] pre- 
sented a toroidal field consisting of concentric rings and an- 
other axisymmetric field, respectively. Spiral patterns based ei- 
ther on the structure of the NE2001 thermal electron density 
model 11611 or on the spiral structure of the Milky Wa y 1711 de- 
duced from HII regions and giant molecular clouds [18] have 
also been proposed to describe this field. The implications on 
the GMF modelling of recent rotation measure maps were re- 



ported in [ 19, 20, 21]. Some of the most recent constraints on 



the regular disk field have been presented in ||22| 1231 |24J . 

The authors of Ref. 1I25I1 presented first a GMF model adding 
a halo contribution made of toroidal and poloidal fields to the 
disk field 11251 12611 . In the present paper, we refer to this model 
as the "PS model". Several new GMF models were proposed 
and confronted with the data in Refs. 11271 l28l 12911 . At present, 
no theoretical GMF model can fit all experimental data IBOlliTll . 

Earlier works discussing the propagation of UHECRs in tur- 



bulent fields include Refs. I32tl33ll . The first studied the generic 
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theoretical properties of multiple image formation and ampli- 
fication in a turbulent magnetic field for point-like UHECR 
sources 13211 . Ref. 13311 discussed the range of values one could 
expect for the relative ratios of UHECR deflections in the reg- 
ular and turbulent GMF. In our present work, we deal with the 
case of iron nuclei at the highest energies emitted by extended 
extragalactic sources considering both the regular and the tur- 
bulent components of the GMF. 

Our paper is structured as follows: Section |2] presents and 
discusses the two numerical methods that we use to generate 
the turbulent GMF. In Section [3] we investigate qualitatively 
the consequences induced by the turbulent GMF on the shape 
of images of single sources and of the supergalactic plane. The 
implications of this field on the (de-) magnification of the fluxes 
from single sources are presented in Section |4] In Section |5] 
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we rule out the possibility that all presently observed UHECR 
events could have been emitted by only one nearby source 
and deflected over the celestial sphere by a strong turbulent- 
dominated GMF. 



2. Modelling of the turbulent Galactic magnetic field 

In this section, we present the two methods we use to model 
the turbulent Galactic magnetic field and recall some of its ba- 
sic properties. A turbulent magnetic field B satisfies (B(r)) = 
and ^B 2 (r)^ = B 2 ms > 0. At the energy we consider throughout 
this work, E — 60 Ee V, the Larmor radius of iron primaries is 
typically much larger than the coherence length L c of the turbu- 
lent field. As a result, the effect of a uniform turbulent field can 
be parametrised in this regime by L c and B lms . Therefore, our 
results for the deflections are largely independent of the exact 
shape of the fluctuation spectrum which is poorly constrained. 
In the following we will consider for definiteness a Kolmogorov 
spectrum. If k denotes the modulus of wave vectors, the power 
spectrum of such a field is Vik) oc k ' 3 and the amplitudes B(k) 
of its Fourier modes follow |B(k)| 2 oc A: -11 '' 3 . 

The modulus of wave vectors k in Fourier space satisfy 
27r/L m ax < k = |k| < 27!-/L min , where L min and L max denote 
the minimal and the maximal scales of variation present in the 
field B(r), respectively. The correlation length L c of the field is 
equal to 



1 1 (^min/^max) 

5 1 — (L m [ R / X max ) 2 / 3 



(1) 



Turbulent magnetic fields can be generated either directly in 
coordinate space as a superposition of a finite number of plane 
waves, or by computing the field on a three-dimensional grid 
using the Fast Fourier Transform (FFT). In the present paper 
we use the first method in Sections[3]and[4]following Ref. [34]. 
The field B(r) is written as the sum of N m modes, 



computed with Eq. (f2]i at any point in space with arbitrary pre- 
cision for the given realization. Its drawback is that computa- 
tions are much slower than with the second method described 
below, because the method requires to compute many trigono- 
metric functions at each step of the particle trajectory. 

The second method consists in precomputing the field on a 
three dimensional grid in real space. The coordinates of its 
wave vectors correspond to positions of vertices on the recip- 
rocal grid in the Fourier space. The resulting field in real space 
is computed with the Fast Fourier Transform (FFT) [35 1 . The 
magnetic field can be quickly interpolated at any point of the 
space. We take an 8 point-interpolation of the field values on 
the eight vertices of the grid cell in which the considered point 
lies. 

We use this method in Section to derive results which re- 
quire the propagation of 3.6 x 10 6 particles. For this work, us- 
ing a three dimensional cubic grid of N = 2 s vertices per side 
is sufficient. It contains 256 3 values for each component of the 
turbulent field. This grid of 2.6 kpc lateral size is periodically 
repeated in real space. 

In the three following sections, we use the same type of 
generic profile for the rms strength of the turbulent field, 
B r ms(r,z)- We take the model introduced in Ref. [36], 



(5) 



B rms (r,z) = B(r)exp 

Zo 



B(r) = ^ A„e„ exp (ik„ ■ r + iB n ) 



(2) 



where r is the galactocentric radius and z the distance to the 
Galactic plane. The parameter zo characterizes the extension of 
the random field in the halo. In the next sections, we will vary 
Zo in the range 0.75 kpc to 10 kpc. The radial profile B(r) is 
defined as 

f Bq exp (|f ) , if r < 3 kpc (bulge) 

B(r) = { / -( r -8.5kpc) \ (6) 

[ B ° eX P( 8.5 kpc J ,lf>>3k P C 

where Bo denotes the value of B rms close to the Sun. For the 
regular field, we use the PS model choosing the parameters as 
described in Ref. 



where s„ = cos(a„)x„ + i sin(a„)y,,, a„ and/3,, represent random 
phases and the unit vectors %, andy,, form an orthogonal basis 
together with the randomly chosen direction of the n-th wave 
vector k„. The amplitude A„ of the n-th plane wave is given by 



a2= BLG(k) 



££i G(k„) 



with 



G(k) = 



4nk 2 Ak 
1 +(fcL c ) 11/3 



(3) 



(4) 



We choose k n between k m i n = 2n/L max and k max - 2n/L m i n 
with equal spacings in logarithmic scale: A log A; = Ak/k = 
const. The number of modes is set to N m = 1000. We verified 
that with this we obtain the theoretically expected deflections 
in the small-angle scattering regime with an error smaller than 
- 2%. Once the modes are dialed, the turbulent field can be 



3. Qualitative consequences of deflections in the turbulent 
GMF for source searches 

We investigate here qualitatively the influence of the turbu- 
lent GMF on the images of single sources and the supergalac- 
tic plane emitting 60EeV iron nuclei. Since the characteristic 
properties of the turbulent GMF are still poorly known, we test 
the dependence of the results on its main parameters. 

The shape of proton and light nuclei source images has been 
discussed in Refs. ll3~7l[36Tl . We have shown in Ref. ifioll how 
images of extragalactic UHECR sources would look like at 
Earth for different models of the regular GMF for a heavy pri- 
mary composition. Contrary to proton sources, one expects that 
some iron sources have several images even at the highest en- 
ergies. Moreover some of them may be strongly magnified or 
demagnified 11121,138 We investigate below the impact 

of a non-zero turbulent GMF contribution on such results. 
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Figure 1: Images of the supergalactic plane (left panel) and of the Fornax galaxy cluster (right panel), respectively considered as a ±10° band in supergalactic 
latitude and a 5° radius disk, emitting 60 EeV iron nuclei deflected in the PS regular GMF model. No turbulent magnetic field was added. Dark blue points represent 
arrival directions of cosmic rays. The grey lines stand for the supergalactic plane and the ±10° circles in supergalactic latitude. The black disk represents Fornax. 
The maps are in Galactic coordinates with the anti-Galactic center in the center. Increasing / from the left to the right. 



We use the model discussed in Section [2] for the profile of 
the rms strength of the turbulent field B lm .s- For the regular 
GMF, we use the PS model as a generic example. Images of 
given UHECR sources can strongly vary from one GMF model 
to another. However, global properties of the whole sky such 
as mean deflection angles or fractions of the sky with given 
(de-)magnification [ loS] are similar in most recent regular GMF 
models. In this respect, the PS model is representative for a 
generic regular GMF model. 

Large scale structure MHD simulations agree that extragalac- 
tic magnetic fields tend to be strongest in the large scale struc- 
ture, around the largest matter concentrations [40 4ll 42 43 



4411 . However, they disagree on some aspects such as the fill- 
ing factor distributions {i.e. the fraction of space filled with 
fields above a given strength, as a function of that strength [45]). 
This leads to substantial differences on the sizes of UHECR 
deflections predicted by different models, ranging from negli- 
gible |40[ 41] to more than ten degrees for proton primaries, 
even at the very highest energies 11421 14511 . Despite this, even 
in the case of strong deflections, UHECR would still prefer- 
entially arrive from directions towards the large scale structure 
(LSS) of the galaxy distribution. Since the fields in the voids 
are very small, UHECR arrival directions outside the Galaxy 
should point back to the supergalactic plane, even if they do 
not point back to their sources. Therefore, if sources are lo- 
cated in the LSS and if extragalactic deflections are sizable, it 
is sufficient in a first approximation to study the image of the 
supergalactic plane itself. 

We show in Figs. QJ[4] the 60 EeV iron nuclei images of 
the whole supergalactic plane (left columns) and of a specific 
galaxy cluster, Fornax (right columns). 

For the images of the supergalactic plane, we assume that the 
region containing most of the local distribution of matter has 
supergalactic latitudes bsc between +10 degrees. Cosmic ray 
arrival directions at Earth are represented in the figures by dark 
blue points, and the three grey lines correspond to the super- 
galactic plane {bsc - 0°) and the bsc - ±10° circles. The 
60 EeV iron image of the supergalactic plane in the PS model 
is split into two. The regions of the supergalactic plane near 
both Galactic poles cross blind regions and therefore do not 
contribute to the image observed at Earth, see Fig. |5](upper left 



panel). This cuts the supergalactic plane (SGP) into two parts: 
The part of the plane located at / ~ 120° - 150° mostly pro- 
duces the image at I ~ 0° - 130°, and the part at / ~ 300° - 330° 
is responsible for the image at / ~ 250° - 350°. At energies 
E > 120 EeV, the two images start to merge. In the recent GMF 
models presented in Refs. 1E7I 291. the 60 EeV iron nuclei im- 



ages of the SGP display several common features. The largest 
part of the image(s) also appear in the region I ~ 0° - 180°, as 
for the PS model, and the largest contributor to this image is the 
part of the SGP which is located at I ~ 120° - 150°. In all tested 
models, the overall geometry of the image in the I ~ 0° - 180° 
region is "circle-like" as for the PS model. The image may 
be enlarged and shifted on the sky by more than 30° from one 
model to another, and its details are model-dependent. 

For the galaxy cluster images, we choose the Fornax cluster 
because its high Galactic latitude (in absolute values) avoids 
complicated images due to the magnetic field in the Galac- 
tic disk. In the PS model, 60 EeV iron images from the two 
nearby clusters with higher latitudes, Virgo and Coma, turn 
out to be invisible at Earth due to their strong demagnifica- 
tion jToll . Fornax is regarded here as a source extending over 
a 5° radius. Even in the case the cluster would contain only one 
source, deflections by magnetic fields inside the cluster lead to 
an extended source J46|. For the galaxy cluster images, we 
also assume that deflections of UHECR nuclei in extragalactic 
magnetic fields are weak enough (as in, e.g. the simulations of 
Refs. 1 4(J 41]) that the UHECRs are not spread out over large 
regions of the LSS. In the figures, the cluster is denoted by a 
black disk. 

To generate these figures, we backtrace 10 4 iron antinuclei 
emitted isotropically from the Earth to the outside of the Galac- 
tic halo. We record those which escape the Galaxy in directions 
in the bsc - -10° to + 10° band for the supergalactic plane 
images, and within 5 degrees from the center of Fornax for the 
Fornax images. 

Figure [T]presents the supergalactic plane and Fornax images, 
deflected in the PS regular GMF only. This figure can be con- 
sidered as the reference to which Figs. |2}]4]should be compared. 

For Figs. |2] [3] and |4] a non-zero turbulent GMF component 
was added to the PS regular component. They show the depen- 
dence of the deflections on the turbulent field rms amplitude Bo, 
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Figure 2: Impact of the turbulent field strength at Earth Bo on the 60EeV iron nuclei images of the supergalactic plane (left column), and of the Fornax galaxy 
cluster (right column). Same GMF regular component as for Fig. ^(PS model), with in addition a non-zero turbulent Galactic magnetic field. It has a Kolmogorov 
spectrum of correlation length L c = 50 pc. Its extension into the halo is zo = 3 kpc. See text for details on the field profile. First row: Bo = 0.25 /jG; Second row: 
B = 1 //G; Third row: B = 4//G. Same key as for Fig.^ 



its correlation length L c and its extension above the Galactic 
plane zo- For clarity, we choose a given "reference" configura- 
tion of the turbulent field parameters {Bq,L c ,zq\, and for each 
figure, only one parameter varies, while the two others are fixed 
to the reference values for which we take Bq - \ yt/G, zo - 3 kpc 
and L c = 50 pc. This configuration corresponds to the second 
row in Fig. [2] For all examples presented in this work, we adopt 
a Kolmogorov spectrum with L m i n = 0.4 L c and L max = 4L C . 

Each plot shown in this section is computed for one given re- 
alisation of the turbulent field. We have checked that the results 
would not be qualitatively different if one takes another reali- 
sation of the turbulent GMF. The only exception to this is the 
lower row in Fig. [3] for the turbulent field with L c = 200 pc and 
Zo = 3 kpc. Since zo is not much larger than L max , cosmic rays 
are more sensitive to the given realisation of the field. 

In Fig. [2] each row corresponds to a different strength of the 
turbulent GMF. From the top to the bottom, the field value at 
Earth Z?o is respectively set to 0.25 fiG, 1 fiG and 4-fiG. The ef- 
fect of this parameter is easily visible both in the Fornax and 
supergalactic plane images. When the turbulent field strength 
is increased, the images are more and more spread out, com- 
pared to those in FigQ] For example, the angular spread of the 



Fornax image at (I ~ 0°,b ~ -70°) increases from negligible 
to nearly comparable to the magnitude of the deflection in the 
regular field itself. At the same time, the "center" of the image 
of Fornax stays at the same place as in the regular field. On the 
two last rows, a few cosmic rays are even deflected to the north- 
ern Galactic hemisphere. For the case Bo — 4 fiG, the structure 
of the supergalactic plane image becomes practically unrecog- 
nizable. For this regular GMF model and the parameters of the 
turbulent component considered here, deflections are however 
still not sufficiently large to spread the image over the whole 
celestial sphere. 

Fig. [3] presents the impact of the correlation length L c . In 
the first row, L c is 4 times smaller than the reference value, 
L c = 12.5 pc. Images are less extended and more compact than 
in the case L c = 50 pc. They nearly resemble to the images in 
Fig. Q] computed without any turbulent field. On the second 
row, L c is set to 200 pc. The spread of the images is larger. 

In Fig. [4] we investigate the impact of the turbulent field ex- 
tension into the Galactic halo. The first and second rows respec- 
tively represent the cases zo = 0.75 kpc and zo = 6 kpc. The fine 
structure of the images is indeed more spread for zo — 6 kpc 
than for zo = 0.75 kpc, but the effect is less remarkable than the 
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Figure 3: Columns and key as in Fig. [2] Impact of the turbulent Galactic magnetic field correlation length L c on the supergalactic plane and Fornax 60EeV iron 
images. Turbulent field with a Kolmogorov spectrum, Bo = 1 t*G and zo = 3kpc. See text for details on the field profile. First row: L c = 12.5 pc; Second row: 
L c = 200 pc. PS model for the regular Galactic magnetic field. 
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Figure 4: Columns and key as in Fig. [2] Impact of the turbulent Galactic magnetic field extension into the halo zo on the supergalactic plane and Fornax 60 EeV 
iron images. Turbulent field with a Kolmogorov spectrum of correlation length L c = 50 pc, and Bq = \ /iG. See text for details on the field profile. First row: 
20 = 0.75 kpc; Second row: zo = 6kpc. PS model for the regular Galactic magnetic field. 



dependence on Bo in Fig. [2] 

The mean deflections in the turbulent field grow with Bq, L c 
and zo- In the limit case of propagation in the ballistic regime, 
the rms deflection S lms for particles of rigidities E/Z propagat- 
ing the distance L in a turbulent field of constant rms strength 
B lms is equal to 



1 ZeB n 



(LL C ) 



1/2 



6.3° x 



Z 60 EeV B, 



26 E ljuG\3kpc/ \50pc 



1/2 



1/2 



(7) 



If both the source and its image are located on the same side 
of the Galactic plane and at high enough Galactic latitudes, 
the distortion of the image follows the simple expectation of 
Eq. (O. The image of Fornax at (I ~ Q°,b ~ -70°) is a good 
illustration of this dependence on the turbulent field parameters 
-see Figs. QJ[4] In this case, iron nuclei from one given source 



with E > 60 Ee V are more or less spread around the arrival di- 
rections that they would have had in the regular field only. Be- 
cause of the linear dependence on fi nns , the size of the spread is 
more sensitive to a relative change of Bo than to a change of L c 
orz - 

Contrary to this simple case, if the source or its image is lo- 
cated at "low" Galactic latitudes (\b\ < 30 - 40°), even rela- 
tively small changes of trajectory due to the turbulent field can 
lead to very large distortions of the image. Some images may 
even appear in completely different locations. Such a behavior 

can be seen for the image of Fornax at (7 ~ 15°, b 15°), see 

Fig. [T] Depending on the parameters of the turbulent magnetic 
field, this image is partly -or completely- shifted in other parts 
of the sky, and new images with Galactic latitudes \b\ < 40° 
appear -see Figs. [2|-|4] Such effects cannot be described by the 
simple analytical formula of Eq. (0. A full study requires the 
numerical computations conducted in this work. 

For the ranges of turbulent GMF parameters considered here, 
we find no noticeable qualitative difference on source images 
between a field with a Kolmogorov spectrum and a field with 
the same modulus for all wave vectors k in Fourier space, i.e. 
£min = ^max = 2L C , as long as the fields have the same correla- 
tion length. 

For such parameters, iron nuclei from a source with E > 
60 EeV are not sufficiently spread to cover the whole sky. How- 
ever, if one takes the extreme case of a very strong field, 
Bq = 10 fiG, extending very far into the halo, zo = lOkpc, 
both the images of single sources and of the supergalactic plane 
are spread over the whole sky. Therefore, we investigate in 
Section [5] if the Pierre Auger data published recently [5] may 
be compatible with only one nearby UHECR source, such as 
Cen A, in case of a very strong and extended turbulent GMF. 

4. Flux (de-) magnification of sources and extension of blind 
regions 

In this section, we investigate the impact of the turbulent 
GMF on the (de-) magnification of individual source fluxes. 
Magnification and demagnification effects are due to magnetic 
lensing of UHECR in the GMF. In Ref. we studied the 
contribution of the regular GMF component. We have shown 
that in the case of iron primaries, source fluxes can be strongly 
magnified or demagnified even for energies as high as 60 EeV. 
In particular, we found large regions of the sky which can be 
considered empty for a given number of backtraced nuclei: For 
instance, regions with an under-density p/{p) < 10 -2 ' can en- 
compass one fifth of the sky. Thus sources located in these parts 
of the sky would not be detectable at Earth by current and next 
generation UHECR experiments. We now study how these pre- 
vious conclusions are affected by the turbulent component of 
the GMF. 

The sky maps in Fig. [5] show the logarithm of the flux am- 
plification factor of extragalactic 60 EeV iron nuclei sources, 
depending on their positions on the sky. We take the PS model 
for the regular GMF. We use the same method as in Ref. II 1 CHI to 
compute M: We isotropically inject 10 5 iron antinuclei at Earth 
and backtrace them to the border of the Galaxy. Their relative 



densities on the celestial sphere outside the Galaxy correspond 
to y\. These densities are computed on 3° radius circles cover- 
ing the whole sky. 

Let us call here "blind regions" the regions of the extragalac- 
tic sky in which log 10 (y[) < -2 for a given energy. At this en- 
ergy, sources located in such regions of the sky cannot be seen 
at Earth by current or next generation UHECR experiments, due 
to the strong demagnification of their fluxes. In general, the 
places, extensions and shapes of blind regions do not change 
abruptly with energy but vary continuously with energy. For 
instance the blind regions in the PS model at 80 EeV are similar 
to those at 60 EeV, though slightly less extended. In Fig. [5] we 
plot blind regions in black. The amplification is below 10" 2 in 
some places. 

The upper left panel in Fig.[5]represents the reference plot in 
which the turbulent field is set to zero. In the three other panels, 
a turbulent component with L c = 50 pc is added to the regular 
field. Its strength and extension in the halo are, respectively, 
equal to (Bq = 1 //G, zo - 3 kpc), (Bq = 4 yuG, zo = 3 kpc) and 
(Bq = 10//G, zo = 10 kpc), for the upper right, lower left and 
lower right panels. 

Figure[5]shows that when the amplitude and the spatial exten- 
sion of the turbulent GMF halo are increased, the extension of 
regions of extreme magnification and demagnification tend to 
shrink. For example, the reduction of the size of the darkest and 
brightest regions from the upper left to the lower left panels is 
clearly visible. For the lower right panel (extremely strong and 
extended turbulent field), regions with -0.5 < log 10 (yi) < 0.5 
encompass =s 80% of the sky. This means that in 80% of the 
sky, 60 EeV iron nuclei sources would have their fluxes neither 
magnified or demagnified by more than a factor ^ 3, In the 
case without any turbulent field, this fraction of the sky with 
moderate flux modification falls to 55%. Blind regions even 
disappear in the lower right panel. 

Figure [6] (left panel) shows four histograms corresponding to 
the four cases shown in Fig. [5] The histograms represent the 
fractions of the sky with given amplifications These his- 
tograms are computed for 10 5 backtraced antinuclei with the 
method used for Fig. 8 of Ref. 10 5 antinuclei are sufficient 
to probe the bin log 10 (y[) < -2. This figure confirms the effect 
of the turbulent Galactic magnetic field. As previously pointed 
out, the distributions tend to peak around the mean value 1 
and have lower variances for stronger and more extended tur- 
bulent components. 

In addition, we can compute the maximum and minimum 
logarithmic amplifications log 10 (y[) on the whole sky. For the 
case without the turbulent field, we have log 10 (y[) 1.15 for 
the maximum amplification. When including turbulent fields 
with the parameters (Bo = 1 A'G, zo = 3 kpc), (Bo = 4/iG, 
Zo = 3 kpc) and (Bo = 10 juG, zo = 10 kpc), we respectively find 
that the maximum values for log 10 (y[) are =s 0.95, 0.63 and 
=! 0.47. For the case with Bo = 10//G and zo = lOkpc, the 

minimum value for log 10 (yi) is 1.6. For the other cases, the 

minimum amplifications are too low to be probed with only 10 5 
backtraced antinuclei. 

The larger the amplitude and extension of the turbulent GMF, 
the weaker the amplitude between the maximum magnification 
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Figure 5: Logarithm of the amplification factor J?( of the flux of extragalactic 60 EeV iron nuclei sources, depending on their position on the sky. Plots are in Galactic 
coordinates, with the Galactic anti-center in the center. We take the PS model for the regular component of the GMF. Regions with logm^) < -2 are in black. 
Upper left panel: No additional turbulent field component is added. The small fraction of the sky with log 10 (J?l) > 1 is in white; Upper right panel: Additional 
turbulent field with Bo = 1 pG, L c = 50 pc and zo = 3 kpc; Lower left panel: Bo = 4/uG, L c = 50 pc and zo = 3 kpc; Lower right panel: Bo = 10//G, L c = 50pc 
and zo = 10 kpc. 
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Figure 6: Histogram of fractions of the sky outside the Galaxy with amplifications J?( within bins of width Alog|o(J?I) = 0.5. Left panel: For all four histograms, 
the GMF regular component is given by the PS model; Middle panel: Sun08 model -"ASS+RING" model of Ref. 1 27] ; Right panel: Sun08-MH model, as defined 
in Ref. 1 10]. In each plot, solid red lines for no additional turbulent field component. Green, blue and magenta lines for the cases of a turbulent field with parameters 
respectively set to (Bo = 1 /jG, zo = 3kpc), (Bo = 4//G, zo = 3kpc) and (Bo = 10/jG, zo = lOkpc). The correlation length L c is 50 pc. The bin of amplifications 
below 10~ 2 corresponds to blind regions. 



and demagnification. This result is in line with our previous 
findings. 

The fraction of the sky covered by blind regions corresponds 
to the left-most bin of the diagram in Fig.|6](left panel). Without 
taking into account the turbulent field, we find that blind regions 
encompass 18% of the sky for the PS model. When including 
turbulent fields with the parameters [Bo = 1 yuG, Zo = 3 kpc), 
(B = 4/j.G, zo = 3 kpc), (B = 10yuG, zo = 10 kpc), blind 
regions respectively represent ^ 16%, - 5% and 0% of the sky. 
Therefore, if the turbulent component is strong and sufficiently 
extended in the z-direction, the fraction of blind regions on 
the sky can easily fall below 15%. However, this would not 
facilitate source searches due to the bigger spread and blurring 
of their images, cf. the previous section. 

We test to which extent these results are generic. We show 



in the two other panels of Fig. [6] the same diagrams for two 
other regular GMF models. In Fig. [6] (middle panel), we take 
the "ASS+RING" model of Reference which we will call 
the "Sun08" model. In Fig. |6](right panel), we use the modified 
version of the Sun08 model introduced in Ref. ifioll . We call it 
the "Sun08-MH (Modified Halo)" model. The four diagrams 
correspond to the same turbulent field strengths and extensions 
in the halo, as in the left panel. One can see that the sizes of 
blind regions without turbulent field are very similar in the three 
models, though they decrease faster with increasing turbulent 
field contributions in the case of the PS model. The tendency is 
however generic: With increasing turbulent field strengths and 
extensions, the sizes of regions of extreme magnification and 
demagnification shrink and eventually disappear. The amplifi- 
cation distributions tend to peak more around the J[ ~ 1 value. 



7 




Figure 7: Left panel: Image of iron nuclei with E > 55 EeV emitted by Cen A, deflected in the PS regular GMF model only, and convolved with the Pierre Auger 
Observatory exposure; Right panel: Same as in the left panel, with an additional strong and extended turbulent Galactic magnetic field. The turbulent magnetic 
field has a Kolmogorov spectrum, with L c = 50pc, Bo = 10/iG and zo = lOkpc. Sky maps in Galactic coordinates, with the anti-Galactic center in the center. 
Cosmic rays are in dark blue and the black disks denote the position of Cen A on the sky. The red circles surround the region within 18° from Cen A. 



For all three regular GMF models, the green lines (Bo — 1 /iG, 
Zo - 3 kpc) are still close to the red lines (no turbulent field), and 
the blue lines (Bq = 4 //G, zo = 3 kpc) show a visible decrease 
of the sizes of the extreme (de-)magnification regions, as well 
as an increase of bins with more moderate M. In all three pan- 
els, blind regions completely disappear for the largest turbulent 
field contribution (Bq = 10/zG, Zo = 10 kpc). All distributions 
are centered on moderate M bins, but their widths depend on the 
regular GMF model. The width for the Sun08 model magenta 
distribution is larger than for those of the two other models. 
This is due to the large strength of this field in the halo -up to 
lOyuG, which is comparable to that of the turbulent field. Re- 
gions with -2 < logj 0?{) < -1.5 and with 0.5 < log 10 (^l) < 1 
indeed disappear when a more reasonable regular field strength 
in the halo is taken, such as for the Sun08-MH model. 



5. Can all UHECR measured at Earth come from a sole 
nearby extragalactic source? 

In this section we test the hypothesis that all UHECR de- 
tected at Earth may be due to one nearby source such as Cen A. 
Cen A has been suggested to be a source of UHECR since a 
longtime 

Since we are sufficiently far from the most prominent objects 
of the nearby Large Scale Structure (LSS) and since Cen A is 
only 3.4 Mpc away and located within a relatively cold region 
of the local LSS, the extragalactic magnetic fields are unlikely 
to cause large deflections between Cen A and the Earth. We will 
neglect them in the following. Therefore, we only have to con- 
sider deflections in the GMF. We plot in Fig Q (left panel) the 
image of Cen A emitting iron nuclei with E > 55 EeV, deflected 
in the PS regular GMF model only. The blue dots correspond 
to the arrival directions of such events, and the black disk de- 
notes Cen A position on the sky. The image is convolved with 
the Pierre Auger Observatory exposure. We use the exposure 
formula of Ref. §5(& with 6 = -35.2° for the Auger declina- 
tion and 6 m = 60° for its maximum zenith angle. We backtrace 
antinuclei from the Earth to outside the Galaxy, with initial di- 
rections following the exposure. The antinuclei energies follow 
the energies of the 69 events published by the Pierre Auger Col- 
laboration [5]. Once the antinuclei leave the Galaxy, we save 
those with escape directions on the sky within 5 degrees from 



the center of Cen A. This corresponds to assuming that Cen A 
is a ~ 5° extended source. We could have used more realistic 
geometries than a disk, but this would not significantly change 
our conclusions. 

Reference ^ discussed the 18° radius overdensity of 
events with E > 55 EeV around Cen A, and Ref. 115 ill found 
that when removing this overdense region, the rest of the sky 
above 55 EeV is still compatible with isotropy. This region is 
surrounded with red circles in Fig Q Refs. 11521 15311 suggested 
that the overdensity may be the heavy nuclei image of Cen A. 
Fig Ul (left panel) shows that in the PS regular GMF model, the 
images of Cen A are deflected beyond the =s 18° radius region 
around Cen A. In all other tested models, heavy nuclei are also 
shifted by at least a few tens of degrees from Cen A. However, 
this result is model-dependent and one cannot rule out the hy- 
pothesis of Refs. J52,[53] yet. Reference ll54ll suggested that 
the overdensity around Cen A may be explained by light nuclei 
emitted by Cen A. The distance to Cen A is small enough to 
avoid problems with photo-disintegration. References 115 U 155 ] 
proposed that this overdensity may be the image of Virgo emit- 
ting heavy nuclei, shifted by the regular GMF. 

Fig. [7] (left panel) shows that iron nuclei emitted by Cen A 
cannot be sufficiently spread over the whole sky in regular GMF 
models to explain all Auger data above 55 EeV. Moreover, we 
have seen in Section [3] that for the ranges of turbulent compo- 
nent parameters used there, ~ 60 EeV iron nuclei would not be 
spread over the whole celestial sphere. Therefore, the only po- 
tential way for the one source hypothesis to be true is that the 
turbulent GMF is very strong and extended. We choose as an 
extreme case Bo = 10//G and zo = 10 kpc. We assume the 
turbulent field has a Kolmogorov spectrum with L c = 50 pc. 
These values for the strength and the extension of this field are 
deliberately larger than the maximum values one can find in the 
literature. For example, one of the largest values proposed for 
zo is zo ~ 7 kpc ||56Tl . 

Blue dots in Fig. [7] (right panel) show the arrival directions 
of iron nuclei with E > 55 EeV emitted by Cen A, when 
this extreme turbulent component is added to the PS model. 
The Cen A image is convolved with the Pierre Auger expo- 
sure. Since deflections are dominated by the turbulent GMF in 
this case, using another reasonable model for the regular GMF 
would not noticeably affect the results presented below. Due to 
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Figure 8: Fractions of cosmic rays within a given angular distance (in degrees) 
to Cen A. Blue distribution for isotropic arrival directions, red one for the 69 
Auger events published in Ref |5], and green one for iron nuclei with E > 
55 EeV emitted by Cen A and deflected in the Galactic magnetic field model 
considered for Fig.[7j(right panel). Both blue and green distributions take into 
account the exposure of the Pierre Auger experiment. 



computing time reasons, we used in this section a turbulent field 
generated via the FFT method, see Section [2]for more details. 
In order to have sufficient statistics for the following study, we 
backtraced 3.6 x 10 6 particles from the Earth. 

The turbulent field starts to be sufficiently strong to spread 
the arrival directions of Cen A nuclei over the whole celestial 
sphere. However, even with such a field, the arrival directions 
of cosmic rays still display a non-negligible dipolar pattern on 
the sky. The density of cosmic ray arrival directions is larger 
in the half-sky centered on Cen A than in the opposite half-sky. 
This tends to be in conflict with the angular distribution found 
from the observational data. 

In Fig. [8] we draw the distributions of the fractions of cosmic 
rays within given angular distances from Cen A, both for the 
published Pierre Auger data with E > 55 EeV [5] (N\ =69 
events) and for N2 — 6300 cosmic rays emitted by Cen A and 
deflected in our extreme GMF model. These two distributions 
are respectively represented by the red and green curves in Fig- 
ure [8] The blue curve on the same panel corresponds to the 
case of isotropic arrival directions of cosmic rays (computed 
with N2 — 10 5 events). Both blue and green lines take into ac- 
count the exposure of the Pierre Auger experiment. The green 
line is always above the blue one. This confirms the presence of 
a non-zero dipolar pattern in the arrival directions of the cosmic 
rays emitted by Cen A. 

We verify with a Kolmogorov-Smimov test i35ll that the hy- 
pothesis that all cosmic rays are emitted by Cen A (green dis- 
tribution) is not compatible with the Pierre Auger data [5] (red 
distribution). This test consists in measuring the maximum dis- 
tance d maK between the two distributions. The significance level 
of the observed value d max , P{d mdx > observed), is directly in- 
ferred from afmax, N\ and N2 [35]. The maximum distance be- 
tween the red and green curves is d max 0.28, at an angular 
distance to Cen A c/cenA - 87°. It corresponds to a probability 
P (aUx > observed) ~ (2 - 4) x 10~ 5 . 



For a weaker or less extended turbulent GMF than the ex- 
treme case considered here, events would be less spread on the 
sky, and the probability P{d m . dx > observed) would be even 
lower. Therefore, the hypothesis that Cen A is the sole UHECR 
accelerator seen at Earth at the highest energies is strongly dis- 
favoured by the Pierre Auger data, even if we consider the most 
favourable case of a heavy composition and a very strong tur- 
bulent GMF component. 



6. Conclusions and perspectives 

In the present paper, we have investigated the consequences 
of a turbulent GMF component on the propagation of UHECRs 
for the case of a heavy composition at the highest energies. We 
have backtraced 60 EeV iron antinuclei in GMF models includ- 
ing both a regular and a turbulent component. For the regular 
component we have mostly used the Prouza and Smida (PS) 
model as a generic example representative for regular GMF 
models. We have tested the model-dependence of our results 
by also considering the Sun08 and Sun08-MH models, as in 
Ref. [10]. We have varied the parameters describing the tur- 
bulent field, accounting thereby for the poor knowledge of its 
properties. 

In Section [3] we have computed the 60 EeV iron images of 
the Fornax galaxy cluster and of the supergalactic plane. We 
have shown qualitatively to which extent the turbulent compo- 
nent may spread these images. We have also discussed the de- 
pendence of such a spread on the turbulent GMF parameters. 

In Section |4j we have shown that the presence of a turbulent 
field tends to reduce the extreme (de-) magnification of individ- 
ual source fluxes. We have called "blind regions" those parts 
of the sky in which the flux of UHECR sources is demagnified 
by more than a factor 100. Current and next generation experi- 
ments will not be able to detect sources located in such regions. 
The size of these regions shrinks when including the effect of a 
turbulent GMF: At 60 EeV, for sufficiently strong turbulent field 
strengths of ^ 4 /iG at Earth and large extensions > 3 kpc into 
the halo, the fraction of blind regions was reduced by more than 
a factor 3 in the case of the PS regular GMF model. 

In the last section, we tested the hypothesis that all UHECR 
detected at Earth above a 55 EeV could be due to a sole 
nearby source, such as Cen A. We found that, even in the most 
favourable case of iron nuclei deflected in a strong and extended 
GMF turbulent component, this scenario is not compatible with 
the present data of the Pierre Auger experiment. 

Finally we note that future radio experiments, such as LO- 
FAR and SKA, will improve our knowledge of the properties 
of the GMF turbulent component, such as its strength and its 
extension in the Galactic halo [57, 58]. 
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